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ABSTRACT 

Very recent observations of the 6 Li isotope in halo stars reveal a 6 Li plateau about 1000 times above the 
predicted BBN abundance. We calculate the evolution of 6 Li versus redshift generated from an initial burst of 
cosmological cosmic rays (CCRs) up to the formation of the Galaxy. We show that the pregalactic production 
of the 6 Li isotope can account for the 6 Li plateau observed in metal poor halo stars without additional over- 
production of 7 Li. The derived relation between the amplitude of the CCR energy spectra and the redshift of 
the initial CCR production puts constraints on the physics and history of the objects, such as pop III stars, 
responsible for these early cosmic rays. Consequently, we consider the evolution of 6 Li in the Galaxy. Since 
6 Li is also produced in Galactic cosmic ray nucleosynthesis, we argue that halo stars with metallicities between 
[Fe/H] = -2 and -1, must be somewhat depleted in 6 Li. 

Subject headings: Cosmology - Cosmic rays - Big Bang Nucleosynthesis - Stars: abundances 
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1. INTRODUCTION 

To account for the origin and evolution of lithium, beryl- 
lium and boron, we rely on our understanding of several very 
different aspects of nucleosynthesis, namely: Big Bang, non 
thermal, stellar nucleosynthesis, all of which must be cor- 
related through cosmic and chemical evolution. These rare 
light nuclei are not generated in the normal course of stel- 
lar nucleosynthesis (except 7 Li in the galactic disk) and are 
in fact destroyed in stellar interiors. This explains the rel- 
atively low abundance of these species. While a significant 
fraction of the observed 7 Li is produced in the Big Bang, 
the Big Bang nucleosynthesis (BBN) production of 6 Li, Be 
and B results in abundances which are orders of magni- 
tude below that observed in halo stars. For example, BBN 
production of 6 Li is dominated by the process D(a,7) 6 Li. 
At the baryon density deduced from observations of the 
anisotropics of the Co smic Microwave Ba ckground (CMB) 
radiation by WM AP JSoergel et al. 1 120031). its BBN value 
is 6 Li/H ~ 10" 14 fThomas. Schramm. Olive. & Fieldslll993l 
Vangi oni-Flam et al. 1 119991) . On the other hand t he BBN 
mean value of the 7 Li abundance is, according to Cvburt 
( 12001 7 Li/H = 4.27+J;^ x 10" 10 , according to lCuoco et aLl 
J2001 7 Li/H = 4.9+};^ x 10" 10 , or according to ICoc et al. I 



(2004) 7 Li/H = 4.15^ x 10" 
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As such, the 7 Li/ 6 Li ratio 

in BBN is about 4 x 10 4 . 

The very low abundances of the 6 Li, 9 Be and 101 'B 
isotopes predicted by BBN theory imply that their most 
plausible production process was the interaction of Galac- 
tic Cosmic rays (GC Rs) with the interstellar m edium 
(for a review see iVangioni-Flam. Casse. & Audouzel l2000). 
Of these isotopes, 6 Li is of particular interest be- 
cause it has only recently been measured in halo stars 
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JSmith Lambert. & Nissen lU993| jHobbs & Thorburn1[ T99l 
19971 ISmith Lambert. & Nissen Ill998l ICavrel et ai*Tll999 
[ Nissen etal I 119991 120001 lAsplund et al. I I200l'r i2004b: 
Aoki et al. 2004) thus offering new constraints on the ver y 
early evolution of light elements ( Steig man et al. I Il993l) . 
Many studi es have fo llowed the evolutio n of 6 Li in our Galaxy 
(see e.g. lFields & Olive Hl999bHVangioni Flam et al. lll999l) . 
Of particular importance in this context is the a+a reaction 
that leads to the synthesis of this isotope (as well as 7 Li) and is 
efficient very early in the evolutionary history of the Galaxy. 

Different scenarios have been discussed to explain the 

abundance of 6 Li in metal-poor halo stars (MPHS). 

ISuzuki & Inoue I (12002) discussed the possibility of cosmic 
rays produced in shocks during the formation of the Galaxy, 
which was consi stent with 6 Li data available at that time. 
Jedamzik (2000) considers the decay of relic particles, dur- 
ing the epoch of the big bang nucleosynthesis, that can yield 
to a la rge primordial abundance of 6 Li. iFields & Prodanovic I 
(2004) have studied in detail the lithium production in con- 
nection to gamma rays, using a formalism similar to ours but 
with a different point of view as fa r as the observational con- 
straints are concerned (see Section l5~2l . 

Until recently, the abundance of 6 Li had been observed in 
only a few MPHS with metallicity [Fe/H] larger than -2.3. 
New values of the ratio 6 Li/ 7 Li have been measured with 
UVES at the VLT-UT2 Kueyen ESO telescope, in halo stars 
with metallicity ranging from -2.7 to -0.5 (see Section |2jl. 
These observations indicate the presence of a plateau in 6 Li/H 
~ 10 _1 1 which suggests a pregalactic origin for the formation 
of 6 Li. 

In this paper, we consider the synthesis of lithium due to the 
interaction of cosmological cosmic rays (CCRs), produced at 
an early e poch, with the intergalactic medium ( Mon tmerle I 
Il977albld and Section [3}- As a + a processes also produce 
7 Li, these models are constrained by the 7 Li plateau observed 
in the same MPHS. This constraint is made more severe by 
the current discrepancy between the BBN predicted value of 
7 Li and the observational abundance. We demonstrate how 
this model can explain the recent observations and constrain 



2 



the history of cosmological structure formation (Section [4} 
and the galactic evolution of 6 Li (Section [5.1l >. We compare 
these results with the expected evolution of 6 Li from GCR nu- 
cleosynthesis. Without the pregalactic production of 6 Li, the 
latter model can not account for the elevated 6 Li abundances 
at very low metallicity. In contrast, models for which O/Fe 
increases at low metallicity are able to produce sufficient 6 Li 
at low metallicity, without pregalactic production. However, 
in this case, the bulk of the 6 Li data seen in higher metallicity 
stars must be argued to be depleted. The same is true for our 
model of CCR nucleosynthesis, but to a lesser extent. We ar- 
gue that 6 Li data in stars with [Fe/H] = -3 to -4 will be required 
to distinguish between these scenari os. O ur predictions will 
be compared to other work in Section l5~2l and our conclusions 
are given in Section|6] 

2. OBSERVATIONAL AND NUCLEAR DATA 

The determination of the 6 Li abundance in MPHS is ex- 
tremely difficult and requires high resolution and high signal 
to noise spectra due to the tiny hyperfine splitting between 
the two lithium isotopes. The line splitting is only seen as 
the narrowly shifted lines are thermally broadened. Though 
the fits to the width of this feature are sensitive to the 7 Li/ 6 Li 
ratio, it is very difficult to obtain accurate measurements of 
the isotopic ratio. 6 Li can only be realistically expected to be 
observed in stars with high surf ace temperatures and at low 
metallicities of [Fe/H] < -1.3. iBrown & Schramm I dl988l) 
determined that only in stars with surface temperatures greater 
than about 6300 K will 6 Li survive in the observable surface 
layers of the star. At metallicities [Fe/H] > -1.3, even higher 
effective temperatures would be required to preserve 6 Li. 

As noted above, the previous sets of data on the 
lithiu m isotope ratio h as been s ignificantly ex panded by 
lAsplund et al. I J2004bl) (see also iLambertl 120041) with the 
observations of 24 MPHS. Previously, only 3 stars with 
metall icity [Fe/H] < —1-3 showed net detections of 6 Li 
JSmith. Lambert. & Nissenlll993[lHobbs & Thorburnl| 19941 
1 19971 ISmith Lambert. & Nissen I [19981 ICavrel et ai"Tll999t 
Nisse net al. 112.00(1) . The observed abundances of Li/H and 
6 Li/H are displayed versus the metallicity, [Fe/H], in Fig.^ 
There are in addition several stars with metallicities in the 
range [Fe/H] = -3 to -0.5 for which only upper limits (not 
shown) to the 6 Li/ 7 Li ratio are available. Note that the 6 Li 
abundance at solar m etallicity is plotted for both the mete- 
oritic value tLodde rs Il2003 |) and the s olar photospheric value 
( Asplund, Grevesse, & Sauval 2004). The latter is derived 
from the photospheric value of Li assuming the solar ratio 
7 Li/ 6 Li =12 and therefore really represents an upper limit to 
the 6 Li photospheric abundance. 

These new data at low metallicity reveal the existence of a 
plateau for 6 Li, whose abundance is about 1000 times higher 
than that predicted by BBN. As such, another production 
mechanism which is capable of producing what appears to be 
an initial enrichment of 6 Li in the intergalactic medium is re- 
quired. Here, we concentrate on the interaction of a particles 
present in CCRs produced at high redshift, with He at rest in 
the IGM, as a potential description of this pregalactic enrich- 
ment process. The abu ndan ces in higher metallicity stars will 
be discussed in Section lBTl 

Note that there is, however, considerable dispersion in the 
data, and as noted above, there are many stars for which there 
was no detectable 6 Li, indicating that depletion may have 



played a role in the observed 6 Li abundance for stars in the 
plateau as well. This is in contrast to the 7 Li plateau which 
shows very little dispersion and for which we exp ect the role 
of depletion to have been minor ( Rvan et al. 2000). 

Galactic production of lithium arises from the interaction 
of GCRs with the ISM via the a + a reaction. This is a pri- 
mary process which yields a logarithmic slope of 1 in the 
lithium abundance versus metallicity ([Fe/H]) relation as seen 
in Fig.^ The amplitude for 6 Li product ion i s constrained by 
the abundances of Be and B (see Section lSTl for an additional 
discussion). As one can see, this model can not explain the 
elevated 6 Li abundances at low metallicity and requires some 
6 Li depletion at higher metallicity (as will all of the models 
discussed h ere). 

Recently, Me rcer et al. I J2001I) have performed new mea- 
surements related to the a + a reaction and provide a new 
fit for the production of 6 Li and 7 Li. The calculated 6 Li 
abund ances from the GCR process using the often applied 
Read & Viola ( 1984) cro ss sections an d the fit at higher en- 
ergy provided by Mercer et al. I (|T001) resulting in slightly 
less (30 - 50 % at solar metallicity) 6 Li are compared in 
Fig.[T] In what follows , we use the most recent cross sections 
of M ercer etal.N200 11) . 
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FIG. 1 .— The evolution of Li/H vs [Fe/H]. Here, the evolution of Li/H 
is modeled by GCR nucleosynthesis alone. Predictions using cross sections 
from Mercer et al. 1 2001) (lower line) are reduced compared to those from 
Read & Viola] 119841) (upper dashed line). The abundances of 6 Li in low 
metall icity stars revea l a plateau. The solar abundance of 6 Li from mete- 
orites i Lodders 2003 ) and the upper limit to the solar photospheric abun- 
dance I Asplund. Grevesse. & Sauval 2004) are also shown. 

3. CCR PRODUCTION OF LITHIUM IN THE IGM : 
FORMALISM 

3.1. On the Existence of Cosmological Cosmic Rays 

The existence and global properties of cosmic rays in 
the Galaxy are often related to supernova explosions and/or 
gamma-ray bursts, in massive stars. Motivated by the 



3 



WMA P resul t s indi cating an early epoch of reionization, 
Daien eet al. I (|2004) have developed models that include an 
early burst of massive stars with several possible mass ranges, 
capable of reionizing the intergalactic medium, while satis- 
fying observational constraints on cosmic chemical evolution 
in pre-galactic struct ures and in the intergalactic medium. In 
particular, Daigne et al. (2004) have demonstrated that the 
presence of massive stars (M ~ 40-100 M Q ) is required at 
high redshift (z> 15-20). This early population of stars (pop 
III) is able to reionize the intergalactic medium and generate 
a prompt initial enrichment (PIE) in metals. It is likely that 
particles will be accelerated within the same process. 

Gamma-ray emission, as well as cosmic ray s, may 
also come from active (Stecker & Salamon 1993; 
Mukheriee & Chiang 1999) and normal (Pavlidou & Fields 
2002) galaxies (see also iLemoine II2002I) . Depending on the 
strengh of the magnetic fields in those structures, cosmic 
rays will be confine d or will propagate into the in t ergalactic 
medium (e.g. Berezinskv, Blasi & Ptuskin 1997; Zweibel 
In addition, recent numerical simulations have shown 
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that the formation of large scale structures leads to accretion 
shocks in the baryonic gas, and thus to particle acceleration 
directly in th e intergalactic medium ( Kan g~&" Jones 1 120021 
iMiniati 1120021 iKeshet et al. 120031 iRvu et al. 120031) . Finally, 
at ultrahigh energies, more exotic sources of cosmic rays have 
also been studied (Bhattachariee, Hill & Schramm H99l 
iSiel et al. Ifl999t) . Clearly, there are several viable mecha- 
nisms for the production of CCRs and just as clearly, there is 
a great deal of uncertainty surrounding their production. 

In this paper, CCRs are assumed to be produced in a single 
burst correl ated to a very early g eneration of pop III stars as 
discussed in Daigne et al. 1 1|2004) at a given redshift z s . Note 
that very little is known about the cosmic ray injection spectra 
at these ene rgies. Here, our form alism is directly derived from 
the work of lMontmerlel dl977al) . hereafter M77. We briefly 
summarize this formalism and note explicitly our differences 
with this model. A power-law distribution in particle energy 
is adopted for the CR injection spectrum, 



<t> a {E) = &\2.5K oip {E+E ){E{E + 2E Q )} 
cm~ 2 s _1 (GeV per nucleon) -1 , 



-( 7 +l)/2 



(1) 

which i s the form expected from sta ndard shock acceleration 
theory (Blandford & Eichierl [T987h . & is a normalization 
factor which is fixed by t he value of the injectio n spectral in- 
dex, chosen to be 7 = 3 (Suzuki & Ino ue 120021) . and by z s . It 
will ultimately be constrained by the observed abundance of 
6 Li in the MPHS (see Section^}. E is the kinetic energy per 
nucleon, Eq = 939 MeV is the nucleon rest mass energy and 
K a p = 0.08 is the abundance by number of 4 He/H. Lithium 
production is sensitive to a's with energy E » 10 MeV/n. 

3.2. Transport function in an expanding universe 

The initial burst of cosmological cosmic rays evolves in the 
framework of an expanding universe with a cosmological con- 
stant. 

If Ni(E,z) is the comoving number density per (GeV/n) of 
a given species at a given time or redshift, and energy, we 
define Nin(E,z) = N(E ,z)/«h(z), the abundance by num- 
ber with respect to the ambient gaseous hydrogen (in units 
of (Gev/n) -1 ). The evolution of Nj.n is defined through the 
transport function 



— + dE {bNi * ) + ^ = Q 'V 



Q is a source function which accounts for different sources 
of particle production while 7b is the lifetime against de- 
struction, b describes the energy losses due to expansion or 
ionization processes ((Gev/n) s ). The energy and time de- 
pendencies can be separated as b(E,z) = —B(E)f(z). We 
can distinguish two cases depending on whether losses are 
dominated by expansion or by ionization. The general form 
for the redshift dependence, whe n expansion dominates is 
f R (z) = (1 +z)~ l \dz/dt\Hq l (e.g. IWick. Dermer & Atovanl 
2004). Other contributions to B or /, do not depend on the 
assumed cosmology and are given explicitly in M77. 

Two important quantities, z*(E,E' ,z) and E s (E,z) are used 
in this formalism. Given a particle (a or lithium) with an 
energy £ at a redshift z, z*(E,E',z) corresponds to the red- 
shift at which this particle had an energy E' . E s (E,z) is 
the initial energy required if this particle was produced at 

the redshift of the burst, z s - In particular, z*(E,E s ,z) = z s - 
The equation that defines z* (Eq. A5, M77) is dz*/dE = 

-\B(E)f{z) |dz/df|] _1 (dz*/dz). M77 gives analytical solu- 
tions for z* when ft\ = 0. When SI a ^ 0, it cannot be solved 
analytically when ionization dominates. Integration of this 
equation shows that z*(E,E s ,z) is the solution of 

r E * dE" 

dz"f(z")(\dz/dt\) z „ = j E - — (3) 

We solve this relation for z* numerically whenever analytical 
solutions are not available. 



3.3. The CCRflux and the lithium abundance 

The evolution of the CCR a energy spectrum is derived, 
using Eq. A8 of M77 and the single burst properties, as 



4> a {E) (3 4> a {E s ) 
< P' ME) 



exp R) 
\b{E,z s 



1 



\dz*/dE'\ E , 
(4) 

where <fr a ii (E 7 z) = $ a {E ,z)/«h(z) is the flux of a's per co- 
moving volume 



$ a p(E,z)=pN a , H (E,z) 



(5) 



and (3 (/?') is the velocity corresponding to energy E (E' s ); £ 
accounts for the destruction term (Eq. A9, M77). 

The abundance by number of lithium (/ = 6 Li or 7 Li) of en- 
ergy E, produced at a given redshift z, is computed from 



dN LH (E,z) 
dt 



= / <J aa ME,E')n He (z)<f> a ,H(E' \z)dE' 



= a l (E)K a M4E,z) [(Gev/n)" 1 s 



1 „-ii 



(6) 



(2) 



where a aa ^i(E ,E') = tj[(E)8{E-E/jA). The cross sections 
used have been discussed in Section|2J Note that this equation 
does not take into account the destruction of lithium in the 
intergalactic medium. We show below that this is a reasonable 
approximation. 

Furthermore, we want to compute the abundance of lithium 
in the gas that is present at the redshift of the formation 
of the Galaxy (see below). We assume that all the lithium 
produced will be thermalized in the protogalaxy before stars 
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form. Thus, the quantity that should be compared to the data 

is, 

[//H](z) = [//H]bbn+/ / 9Nl ' H ^ ,z!) dE\dt/dz'\6z'. 

(7) 

where [//H]bbn is the primordial abundance predicted by 
BBN. 

The redshift evolution of 6 Li/H and 7 Li/H are the main re- 
sults that will be compared to observations, and used to con- 
strain the CCR proton energy density, 

4(z)=/ [<S>a{E p ,z)/K ap }E p //3dE p . (8) 

where E cut = 10 MeV corresponds to the a energy cut-off 
(MeV/n) for the a + a — > Li reaction. 

Finally, we comment on the subsequent destruction of 
Lithium. The differential rate of destruction (by protons) is 
equal to cjd (E ) «h (z) Ni h (E , z)(3{E ) and is proportional to the 
Lithium abundance. The cross section ctd, decreases rapidly 
with energy below 10 MeV/n (see Fig. 2 of M77). Assuming 
a constant energy of 10 MeV/n we can derive an upper limit to 
the destruction process. We find that taking destruction into 
account increases the final proton energy density only up to 
7% for z s = 100 and has virtually no effect for z s < 50. 

3.4. Updated quantities 

Since 1977, the cosmological parameters, the 6 Li and 
7 Li abundances predicted by BBN and observed in MPHS 
have changed considerably. They have been updated 
here. Unless oth erwise noted, we use the standard 
ACDM cosmology (Sne rgel etal") 120031) for which H = 
71 kms^Mpc" 1 , Q m = 0.27, Q A = 0.73 with f2 b h 2 = 
0.0224. The Hubble constant is defined as H(z) = 

H (n in {l+z) 3 + {l-n m -n A )(l+z) 2 + n A ) - 5 valid in a 
post-radiation dominated universe and we assume an equa- 
tion of state parameter for the dark energy, w = — 1. Finally, 
dz/dt = -(l+z)#(z). 

We use the recently calculated BBN abundances of lithium 
by ICoc et al. I ((2004) and the observed 6 Li and 7 Li abun- 
dances in MPHS. From Fig. we define the MPHS abun- 
dance for 6 Li as the value corresponding to the plateau, [ 6 Li] p 
= log( 6 Li/H) + 12 = 0.8. 

4. CONSTRAINTS ON CCR PRODUCTION 

The process described above occurs in the IGM, and modi- 
fies the abundance pattern of the medium that will later form 
the Galaxy. Observations of MPHS trace the evolution of the 
gas in the halo of the Galaxy at an early epoch at low metal- 
licity. The observed abundance for the lowest metallicity is 
then assumed to be a pregalactic abundance, i.e. the predicted 
abundance at the redshift of the formation of the Galaxy. This 
may be justified by the presence of the 6 Li plateau. We will 
assume that the pea k of the formation of the structures oc- 
curs at Za^ ~ 3 (e.g. lFontana et al. ll999llJuneau et al. "12001 
Hopkins 2004). Therefore, the abundance observed for the 
lowest metallicity stars must correspond to the abundance in 
the intergalactic medium at z = z ga i- We next define the pro- 
cedure used to constrain the CCR burst parameters from the 
lithium observations. 

4.1. Procedure 



We begin by working within the context of the standard 
framework described in Section|5] that is a ACDM + WMAP 
cosmology. The shape of the CCR spectrum is given by Eq.^ 
with 7 = 3. Then, the evolution of the lithium abundance with 
redshift in our model is uniquely specified by the normaliza- 
tion constant & and the redshift of the CCR burst, z s . 

Our CCR spectrum is constrained by (?) the Spite plateau 
dSnite & Si)itell982l) fo r 7 Li and (z'z) the hint for a 6 Li plateau 
( Asplund et al. 2004b). The Spite plateau should correspond 
to the primordial value of the Li abundance. However, the 
observed 7 Li abundance is a factor of 2-3 lower than the cal- 
culated one (based on the WMAP baryon density). As a result 
of this discrepancy, the 7 Li plateau acts as a strong constraint 
in our model, since it forbids us to produce a non-negligible 
amount of Li. This constraint is weakene d if the observa- 
tional value of 7 Li were higher (see e.g. Mel endez & R amirez 
2004) as the GCR component of 7 Li would become more dif- 
ficult to obser ve as the ratio of GCR to BBN produ ced 7 Li 
is diminished ( Fields. Olive & Vangioni-Flam 2004). Nev- 
ertheless, our model must produce a small quantity of this 
isotope compared to the BBN abundance. The abundance of 
6 Li observed at very low metallicity is assumed to trace the 
abundance in the intergalactic medium before the formation 
of the Galaxy. As mentioned above, our model must be able 
to reproduce this abundance at z = z ga i. 

The constraints on the 6 Li abundance from the calculated 
BBN abundance at z ~ °° and from its observed 'pregalac- 
tic' value at z = z ga i specify a unique amplitude for the CCR 
energy spectrum, jF, for a given z s - Therefore, the normaliza- 
tion constant, J?, is determined by the choice of parameters 
(f^A, 7j z s ), under the constraint given by the initial and 
final 6 Li abundances. The initial 7 Li abundance is fixed by 
the BBN. Then, for each set of parameters, we check that the 
model does not produce too large of an additional pregalactic 
component of 7 Li. 

4.2. The 6 Li plateau and the evolution of 6 Li versus redshift 

The evolution of the abundance of 6 Li with redshift is shown 
in Fig.|2]for three values of z s =10, 30 and 100. In each case 
a 6 Li plateau is produced. By construction, the abundance of 
6 Li at z = Zgai is fixed. In addition, the rate of production of 
6 Li decreases rapidly soon after the initial burst. This was 
noted in M77 and corresponds to the dilution of the CCR flux 
with the expansion of the Universe. Unless the burst occurs 
just prior to the formation of the Galaxy (Section l5~TT l. the 6 Li 
abundance is almost constant for z < z ga i. 

Since the production rates of 6 Li and 7 Li are similar, the ad- 
ditional production of 7 Li due to CCRs (~ 10 -11 ) is negligible 
compared to the BBN primordial values (~ a few xl0~ 10 ). 
Note that the predicted abundance in the Spite plateau is in- 
creased by only 6%, 8% or 10% for z s =100, 30 or 10 respec- 
tively. The ratio 7 Li/ 6 Li follows the same trend as 6 Li (upper 
panel of Fig. [2} and reaches a final value of about 60. 

4.3. Influence of the different parameters on the CCR 
amplitude 

We next investigate the influence of the parameters of the 
model (Oa, f2 m ; 7, z 4 ) on the required amplitude of the CCR 
flux (#■) using the same constraint on the 'pregalactic' value 
[ 6 Li] p = 0.8 at z = z ga i. Results are given in Tabled Note 
that, as mentioned above, the evolution of the lithium produc- 
tion is dominated by the dilution of the CCR flux that roughly 
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FIG. 2. — Redshift evolution of the 6 Li and 7 Li abundances in the inter- 
galactic medium (lower panel), and of the ratio of 7 Li/ 6 Li abundance (upper 
panel). The redshift of the initial CCR burst z s is chosen to be 10, 30 and 
100 and is represented by the dotted, dashed and solid lines respectively. The 
shap e of the C CR energy spectrum is fixed by 7 = 3 and a ACDM cosmol- 
ogy 1 Spergel et al. 2003) is assumed. The initial abundances of the lithium 
isotopes are fixed according to BBN calculations (crosses and horizontal dot- 
dashed lines) while the abundance of s Li is chosen to be 10~" 2 at the redshift 
of the formation of the Galaxy, Z e-A = 3 (circle). These choices fix the ampli- 
tude of the CCR flux (see Section l4*Tl . We find that the primordial abundance 
of 7 Li is increased by less than 10% from z s to z ga i. 



TABLE 1 

CCR AMPLITUDE VERSUS INPUT PARAMETERS 
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-4.1 


-9.4 


(0.0,1.0) 
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-4.0 


-9.3 


(0.7,0.3) 


3.0 


-4.4 
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-3.7 


-10.6 


(0.0,1.0) 


3.0 


-4.1 


-11.4 




2.0 


-3.4 


-10.3 


(0.7,0.3) 


3.0 


-3.5 


-12.2 




2.0 


-2.7 


-10.9 


(0.0,1.0) 


3.0 


-3.0 


-11.8 




2.0 


-2.3 


-10.5 



The redshift of the CCR burst (z s ), the cosmology, and the shape of the energy 
spectrum (7) are varied. The amplitude of the CCR spectrum is deter- 
mined in order to reproduce the observed 'pregalactic' (z = z ga i) abundance 
[ 6 Li],, = 0.8 (Fig.|2j. g p is the total initial energy density of protons in the 
CCRs(Eq.|8j. 



follows a (1 +zf law. Thus, the cosmological parameters and 
the shape of the energy spectrum have very little influence on 
the shape of the curves in Fig. [2] Only the amplitude of the 
initial CCR flux varies. 

As one can see from the Table, our results are very sensitive 
to the redshift of the burst, z s - At high energy, one can show 

that & a ji(E,z) oc 4> a {E) (i+zi") ■ Thus, the required en- 
ergy density at z = z s in the CCR is roughly proportionnal to 
(1 +Zs) 1 ' 5 . The CCR normalization also depends on the shape 
of the energy spectrum, 7. A steeper spectrum (higher 7) fa- 
vors the low energy part of the spectrum, where the lithium 
production peaks. Thus, for a fixed amplitude jF, the abun- 
dance of 6 Li will be higher for 7 = 3 than for 7 = 2. Con- 
versely, for a fixed abundance, J^" must be lower for 7 = 3. 
Finally, there is little dependence on the cosmological param- 
eters, especially for large values of z s . 

We have used the observed plateau of 6 Li to set the amount 
of pregalactic 6 Li production. Then, assuming a given epoch 
for the formation of CCRs, the amplitude of the energy spec- 
trum, J?, and the energy density are fixed (Table [Q. The 



overall range of the proton energy density, S p at z = z s is 



12 - 2 , for WMAP concordance model when 7 = 3 
jq-12.2 more generally. Yet, those CRs may 



10" 10 - 2 to 10 
or 10" 9 - 3 to 

also play a role in heating and ionizing the IGM at high red- 
shift. In fact, when z s = 10, the energy density of 6.3 x 10~ 13 
ergs/cm -3 is marginally consistent with the resulting temper- 
ature of the IGM today. At higher z s , this constraint is far 
less important as the resulting IGM temperature scales as 
S p l (1 +z s ) 4 and since <§ p increases slower than (1 +z s ) 4 - It is 
interesting to note that CCRs were predicted to heat the IGM 
and thus avoid the problem of overcooling in the IGM gas 
(Blanchard. Valls-Gabaud. & Mamon 1992). Furthermore, 
Nath & Biermann ( 1993) have put other constraints on the 
total luminosity of the CRs from the Gunn-Peterson optical 
depth at z = 4.2, the Compton y-parameter and metal enrich- 
ment. Assuming the production of CCRs from galaxies at 
z = 10, they obtained an initial luminosity of about 1 .6 x 10~ 27 
h ergs/cm 3 /s. Alternatively, from the amount of metals ejected 
by SN, they place an upper limit on the cosmic ray energy 
density of 10~ 14 ergs/cm at z = and solar metallicity. At 
z = 10, the metallicity of the local ISM corresponding to the 
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site of the CCR production is about .01 solar. Note that this 
is much larger than the resulting IGM metallicity where 6 Li 
production occurs. At z = 10 and at a metallicity of 0.01 solar, 
their limit is effectively relaxed to S p < 10~ 12 erg/cm 3 . 

Thus, we see that the CCR production of 6 Li is capable 
of explaining the large abundance of 6 Li with negligible pro- 
duction of Li (both relative to the BBN value). Sufficient 
6 Li production is achieved by adjusting the flux of CCRs and 
depends primarily on the assumed redshift of the initial burst, 
z s - Subsequently, the 6 Li abundance remains roughly constant 
until additional 6 Li is produced in the Galaxy through CGRs 
as we discuss in the next section. 

5. DISCUSSION 
5.1. Galactic evolution of 6 Li 

The above model for the CCR production can be thought of 
as a form of prompt initial enrichment if our Galaxy is formed 
hierarchically from previously evolved structures. CCRs pro- 
duce 6 Li and a small amount (about a few % of the BBN 
value) of 7 Li. Subsequently, the abundances of these element 
isotopes as well as all other element abundances are controlled 
by galactic chemical evolution. At this point, we assume only 
that the initial abundances in the gas in the Galaxy correspond 
to those in the IGM at z = z ga i. 

It is certain that 6 Li will be produced in the ISM through 
GCR nucleosynthesis (described briefly below) as this is the 
primary mechanism for the production of 9 Be and 10 B. These 
isotopes will not have been produced in any significant quan- 
tities as the IGM was initially devoid of C, N, and O needed 
for spallation processes. In contrast, the presence of primor- 
dial He allows for the CCR production of Li. 

Cosmic rays produced in the early Galaxy will invariably 
interac t with the existing ISM. In stan dard GCR nucleosyn- 
thesis (Reeves, Fowler, & Hovle 1970) LiBeB nuclei are pro- 
duced by spallation when protons and as in the cosmic rays 
impinge on ISM C, N, or O. LiBeB is also produced when 
CNO in the cosmic rays are spalled by ISM protons and 
as. As such, spallation requires heavy elements ('metals') 
to be present in either the cosmic rays or the ISM. In addi- 
tion, a + a fusion reactions between cosmic rays and the ISM 
lead to the production of the lithium isotopes. Indeed, these 
were precisely the types of process considered above in our 
model of CCR nucleosynthesis. Note th at 7 Li and n B also re- 
ceive contributions from the is- process (Wooslev HHJ1990; 
lOlive et al. 1119941 l^ngioni-Flam et al. lfl996l) . but these will 
not be important for our present discussion. 

As in the case of CCR nucleosynthesis, a fusion in GCR 
nucleosynthesis is a primary process in contrast to the pro- 
duction of Be and B in standard models. The spallation 
of ISM CNO is a secondary process so these abundances 
scale as the square of a metallicity tracer such as O or 
Fe if [O/Fe] is constant at low values of [Fe/H]. Moti- 
vated by the observational fact that the log of the Be and 
B abundances appear to scale linearly with [Fe/H] it has 
been proposed that the bulk of cosmic rays are not ac- 
celerated in the general ISM , but rather in the metal-rich 
interiors of superbubbles ij C. asse. L ehouca. & Vangion i-FlarrJ 
1995; Parizot & Drurv 1999), specifically at low metallicity. 
Because the superbubble composition is enriched in metals, 
any cosmic rays which are accelerated in superbubble interi- 
ors would have a composition which is both metal-rich and 
time-independent. The low-energy component of the hard en- 
ergy spectra associated with superbubbles results in the pri- 



mary production of Be and B. We refer to this process as LEC. 
In general both the LEC and standard GCR nucleosynthesis 
are responsible for the observed LiBeB abundances. 

In Fig[0 we display the GCR production of 6 Li in the ab- 
sence of the prompt initial enrichment produced by CCR nu- 
cleosynthesis. Here, we have normalized the flux of galactic 
cosmic rays so as to correctly reproduce the solar value of 
Be/H. The overall flux is the only parameter available in GCR 
nucleosynthesis, and as a consequence the abundances of 10 B 
and 6 Li are predictions of the model (recall that the n B and 
'Li receive an additional contribution from the ^-process). 
As exp ected the logarith mic slope of [Lil versus [F e/H] is 1 
( Fiel ds & Olive ll999btlVangioni-Flam et al. 119991) . 

In Fig.|3]we show the evolution of 6 Li vs [Fe/H] when both 
CCR and LEC processes are included. As one can see in Fig- 
ure 3, without the initial enrichment of 6 Li due to CCRs, the 
evolution of 6 Li resembles that of standard GCRs. In this 
case, the 6 Li abundance begins at very low values and rises 
with a slope of unity until late times. This model alone can 
not explain the observational data. At low [Fe/H] (< -2), the 
observed 6 Li abundance is too high to be accounted for by 
standard CGR + LEC nucleosynthesis. 

In addition, to explain the data at higher [Fe/H] (> -2), 
one must argue that depletion has lowered the abundance of 
6 Li. This is perhaps reasonable as the depth of the convec- 
tion zone is increased at higher metallicity for a fixed surface 
temperature. We note that many of the stars observed only 
reveal upper limits to the 6 Li abundance. That is, in roughly 
15 examples of stars with similar temperatures and metallic- 
ities as those shown, no 6 Li was detected. The lack of 6 Li in 
some stars, coupled with the dispersion seen in the data may 
also indicate that some depletion of 6 Li has occurred in some 
of these stars. Indeed, the difference between the solar pho- 
tospheric and meteoritic values corresponds to a destruction 
of 6 Li of at least a factor of about 200. In this model, we 
would argue that the destruction of 6 Li is negligible at [Fe/H] 
< -2 where the calculation from galactic processes cross the 
plateau. 

We also show in Fig.|3]the evolution of 6 Li when the prompt 
enrichment due to CCRs is included. In this case, the data at 
low [Fe/H] is nicely modeled but depletion is still required to 
explain the data at higher metallicity. At present, the evidence 
for the plateau hinges on the abundances in only a few stars at 
low metallicity. However the two models shown in Fig.|3]can 
be distinguished by future observations of 6 Li at metallicities 
[Fe/H] < -2.7 which would establish the role of CCR nucle- 
osynthesis as a mechanism for the early production of 6 Li. 
As mentionned above, a prompt initial enrichment (PIE) in 
heavy elements is also expected in the intergalactic medium, 
especially within this Pop III stars scenario. However, the 
initial mass fraction of ir on may be of the order of X(Fe) = 
10~ 7 dDaigne et al. l l2004) and thus do not modify the curves 
in Fig. H anc l Fig.|3] for [Fe/H] > -4. Thus, models with a 6 Li 
plateau are not affected by the iron PIE generated by pop III 
stars. 

The GCR nucleosynthesis models described above were 
based on the assumption that [O/Fe] was constant at low 
[Fe/H]. That is, we can use the iron abundance to trace 
the evolution of the LiBeB elements. However some 
data show that [O/Fe] increases with d e creasing r Fe/H| 
jIsr^lmri J _Gajcja2L6pe^ ^ Rebolol 119981 iBoesgaard et alJ 
1999; Israelian et al. 2001). As a consequence, the evolution 
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FIG. 3. — As in Fig.Q the evolution of 6 Li/H vs [Fe/H]. In this case the 
evolution of 6 Li/H is modeled by GCR nucleosynthesis with the inclusion of 
the LEC (lower line). Also shown is the case where a prompt initial enrich- 
ment is produced by CCRs in the intergalactic medium (upper line). This 
model can explain the observed 6 Li plateau at low metallicity. The potential 
depletion of 6 Li is indicated by the dashed arrows. 



FIG. 4. — The evolution of 6 Li/H vs [Fe/H]. Here, [O/Fe] is assumed to 
scale with [Fe/H]. We use the maximum slope of the relation that does not 
require a pre-galactic process (see text for details). However, the amount of 
destruction required is much higher, and must be present at metallicity as low 
as -3. 



of Be and B may appear to be primary with re spect to [ Fe/H], 
but in fact is secondary with respect to [O/H] ( Fields & Olive 
1999a). In reality, the data show that with the respect to [O/H], 
BeB have admixtures of primary and secondary components. 
T hat is, the slope for log(BeB/H) vs. [O/H] is between 1 and 
2 JFields etamOOOUKingttOOll). 

In Fig. |3 we show the resulting evolution of 6 Li when the 
slope of [O/Fe] vs. [Fe/H] is taken to be -0.45 (cf . Fields etaD 
2000; [Fields. Olive & Vangioni-Flaml (2004). For a slope 
larger than -0.45, a prompt initial enrichment is not required 
since the galactic production of 6 Li would already exceed 
the abundance observed in all stars. However, within such a 
model, the destruction rate of 6 Li must be non-negligible for 
metallicities as low as [Fe/H]~ -3. In Fig.^ an d Fig-El the 
slope was chosen to be corresponding to constant [O/Fe]. 

Before concluding this part of the discussion, we note that 
the dispersion in the data may be due to irregular production 
rather than depletion. Due to the dilution of the CCR flux, 
the production of lithium saturates soon after the initial burst 
(Fig. at z = z s - However, if this burst does not occur at a 
redshift much larger than z ga i, the abundance of lithium can 
still increase. Fig.[5]shows the expected variation of the 6 Li 
and 7 Li abundances from z = to 3. If the IGM can pollute 
the Galaxy (e. g. through the merging of satellites, see e.g. 
iNavarro 1 12004) at z < z ga i, the abundance pattern in stars that 
form later may reflect this dispersion. For our choice of z ga i = 
3, the late production of Li may constrain the redshift of the 
CCR burst to be greater than about 5-6. 




4 6 8 10 

Zs 

FIG. 5. — Predicted dispersion in the 6 Li and 7 Li abundances from z = to 
3 as a function of z s ■ Note that by construction, Zs must be greater than the 
redshift of formation of the Galaxy, Zgai = 3 (vertical line). 



5.2. Comparison with previous work 

ISuzuki & Inoue ] il200l consider a model where cosmic 
rays generated by structure formation, during the process of 
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Galaxy formation, produce 6 Li in the course of the evolution 
of the Galaxy. In their model, much of the 6 Li is produced 
early and therefore they also predict a 6 Li plateau which ex- 
tends down to at most [Fe/H] s» -3. The characteristics of the 
plateau depend on the history of structure formation. If this 
process oc curs early enough in th e formation of the Galaxy, 
model I of Suzuki & Inoue (2002) is consistent with the new 
observations (see their Fig. 1). Once again, observations of 
6 Li at metallicities between -3 and -4 can distinguish between 
this model and the one we have presented here for which 
the plateau is predicted to extend to much lower m etallicities. 
Furthermore, within the model of Suzuki & Inoue ( 2002), the 
exact evolution of the 6 Li abundance can be linked to the mean 
azimuthal rotation velocity of MPHSs. Consequently, they 
predict larger dispersion among the observed 6 Li abundances, 
together with a correlation between the 6 Li abundance and the 
rotation velocity. This could also help distinguish between 
galactic a+a GCR production from the early a+aCCR pro- 
duction conside red in this paper. 

The work by Fields & Prod anovicl ll2004l) uses a similar 
formalism to that described above, though redshift evolution 
is not formally taken into account. However, their main fo- 
cus is on the lithium-gamma ray connection in relation to 
the solar 6 Li abundance. Under this assumption, they claim 
that, if CCR interactions account for all of the 6 Li produc- 
tion, it will also account for all of the observed extragalac- 
tic gamma-ray background (EGRB). In this paper, we claim 
that CCRs must produce a pregalactic 6 Li abundance, that 
is about 10 times smaller than the solar abundance since 
most of the 6 Li in stars at solar metallicity is produced dur- 
ing galactic evolution (Fig. [Q. Hence, we would argue that 
it should produce only 10% of the total EGRB (from their 
Eq. 13), which is consistent with th e oretical p redictio ns (e.g. 
|Berezinsky 1 Blasi & Ptuskin"| |1997| : [Co lafrances co & Blasil 
ll998UMiniati 120021) . iProdanovic & Fields I ( 12004 argue that 
the observation of Li in high-velocity clouds may help estab- 
lish the necessity of an e arly source of Li. 

Finally, we note that iJedamzik I (120001 12004a b) consid- 
ers the very early produ ction of 6 Li during BBN from decay 
( IJedamzik 1120001 l2004d) or annihilation (Jedamz ikll2004albl) 
of relic particles. Naturally, this model will also predict the 
existence of an elevated plateau. 



iCvburt. Fields. & Oliyd l200l ICoc et al. I l200l iLambertl 
120041 iRvan & Elliot 112004 . On the other hand, the hint for 
a plateau in 6 Li at very low metallicity, and at a higher abun- 
dance than predicted in standard BBN (by a factor of 1000), 
requires an additional process that produces 6 Li in a pregalac- 
tic phase. The process studied in this paper involves the inter- 
action of a particles present in early cosmological cosmic rays 
with primordial Helium present in the intergalactic medium. 
We have shown that it is possible to produce sufficient quan- 
tities of 6 Li, without the additionnal over-production of 7 Li. 

The early production of 6 Li will be present in the gas that 
forms the Galaxy at z — z ga i and provides a simple explanation 
for the existence of the observed 6 Li plateau in MPHS. The 
level of the 6 Li plateau may provide a strong constraint on 
the z s -j£" plane for the initial burst of CCRs, and hence on its 
total energy. However, the existence of this plateau needs to 
be confirmed with additional observations of 6 Li in stars with 
metallicities lower than -3. If the 6 Li plateau persists down 
to lower metallicities, it could confirm the predictions of this 
model and distinguish the physical processes occurring during 
Galaxy formation. 

In a forthcoming paper, we will go further than the sin- 
gle burst approximation. The CCR production could be re- 
lated to the formatio n and chemical evolution of pop III stars 
(Daig ne et al. 12 004). The influence of this process in the pro- 
duction of other elements, such as Be, B and D, will also be 
studied. 
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6. CONCLUSION 

The existence of the Spite plateau for 7 Li indicates that 
low metallicity halo stars are representative of the primordial 
BBN abundance, although the discrep ancy with predictions 
based on WMAP results is still an issue (Aspl undetal 119991 
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